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We study adiabatic pumping of electrons through quantum dots attached to ferromagnetic leads. 
Hereby we make use of a real-time diagrammatic technique in the adiabatic limit that takes into 
account strong Coulomb interaction in the dot. We analyze the degree of spin polarization of 
electrons pumped from a ferromagnet through the dot to a nonmagnetic lead (N-dot-F) as well as 
the dependence of the pumped charge on the relative leads' magnetization orientations for a spin- 
valve (F-dot-F) structure. For the former case, we find that, depending on the relative coupling 
strength to the leads, spin and charge can, on average, be pumped in opposite directions. For the 
latter case, we find an angular dependence of the pumped charge, that becomes more and more 
anharmonic for large spin polarization in the leads. 

PACS numbers: 72.25.Mk, 73.23.Hk, 85.75.-d 
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I. INTRODUCTION 

Charge and spin transport through a nanoscale con- 
ductor can be obtained, in the absence of a transport 
voltage, by periodically varying in time some of its pa- 
rameters. If the time dependence of the system is slow 
compared to its characteristic response time, we refer to 
this transport mechanism as adiabatic pumping. This 
particular regime allows us to study the properties of a 
system being slightly out of equilibrium due to an explicit 
time-dependence of its parameters. Numerous works 
have studied mesoscopic pumps both theoreticall}*ii^i^'^ 
as well as experimentallyj^iii^iSiiS The established frame- 
work to calculate the pumped charge through a meso- 
scopic scatterer is based on the dynamical scattering 
approachflili This approach can be applied when the 
Coulomb interaction can be neglected or treated within 
the Hartree approximation. Recently, the interest in in- 
cluding the effects of Coulomb interaction beyond the 
Hartree level to the problem of adiabatic pumping has 
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Spin-dependent transport through nanostructures has 
recently attracted a lot of interest. A model example is a 
quantum-dot spin valve, which consists of an interacting 
(single-level) quantum dot attached to two ferromagnetic 
leads (F-dot-F), see Fig. [TJ The leads have, in general, 
non-coUinear magnetization directions and different po- 
larization strengths. Transport through a quantum-dot 
spin valve with noncoUinear leads has been studied ex- 
tensively in the dc limit .SAiS^. In magnetic multilayers the 
tunneling current depends chiefly on the relative orienta- 
tion of the magnetization of the ferromagnetic layersi^^iH 
The situation is more complex in a quantum-dot spin 
valve due to the interplay of the lead magnetization. 
Coulomb interaction, non-equilibrium spin accumulation, 
and quantum fluctuations. In particular, a finite spin ac- 
cumulation is generated on the dot, which plays an im- 
portant role in determining charge transport. 





FIG. 1: Schematic setup of N-dot-F or F-dot-F setup. The 
magnetization directions and the polarization strengths of the 
left and right lead can in general differ from each other (left). 
Sketch of the coordinate system and the magnetization direc- 
tions of the leads (right). 



In the present paper we combine the ideas of adiabatic 
pumping and spin-dependent transport through interact- 
ing nanostructures. We consider two scenarios. First, we 
focus on the situation when only one of the two leads 
is ferromagnetic (N-dot-F), for which we study the spin 
pumped into the non-magnetic lead. Spin pumping in 
systems where the spin degeneracy is lifted by a magnetic 
field has been the subject of several studies J^ii^ii^i^^!^ 
Furthermore spin pumping by means of electrical gat- 
ing only was predicted in a system with Rashba spin- 
orbit coupling. Several aspects of a noninteracting s pin 
pump based on ferromagnets were studied in Ref. l28l . 
Spin pumping through an interface between a ferromag- 
net and a non-magnetic metal has been investigated in 
Ref. [2^ where the pumping cycle is realized exploiting 
the precession of the magnetization of the ferromagnet. 
In our setup we are interested in spin pumping obtained 
by varying periodically the properties of the scattering 
region, such as the dot level position and the tunneling 
strength to the left and the right lead, but leaving the 
lead properties, as e.g. their magnetizations, constant in 
time. A particular intriguing result for the system under 
consideration is that, depending on the relative coupling 
strength between dot and leads, spin and charge can on 
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average be transported in opposite directions. 

Second, we consider the case that both leads are spin 
polarized (F-dot-F). We study the influence of the spin 
polarizations of the leads on the pumped charge and on 
the average spin accumulated on the quantum dot dur- 
ing a pumping cycle. As a result, we find that also the 
pumped charge displays the spin-valve effect, i.e., a de- 
pendence on the relative angle between magnetization 
direction of the leads. For stronger spin polarization of 
the leads, the pumped charge becomes a more and more 
anharmonic function of the relative angle. 

In order to calculate the charge and the spin pumped 
through the dot we use real-time diagrammatic technique 
in the adiabatic limi t^^i'^° and perform a rigorous pertur- 
bation expansion in the tunnel coupling to the leads. We 
consider the system in the regime of weak coupling, con- 
sidering only first-order processes in the tunnel-coupling 
strengths. 



II. MODEL AND FORMALISM 

A. Hamiltonian 

We consider a single-level quantum dot contacted by 
tunnel barriers to two ferromagnetic leads with different 
spin polarization axes, as shown in Fig. [T] For finite 
spin polarization in both leads (F-dot-F), the system is 
called a quantum-dot spin valve. The limit of one normal 
and one ferromagnetic lead (N-dot-F) is included by set- 
ting the spin polarization of one lead to zero. The total 
Hamiltonian of the system can be written as 



H = H, 



dot 



E 

Q=L,R 



lcad,( 



tunnel 



(1) 



electron in lead a can take the values cr = ±, where -I- 
refcrs to the majority spin and — to the minority spin of 
this lead. We choose a coordinate system with the three 
basis vectors e^,. By and e^, pointing along riL-l-nR, nL~ 
fiR and Ar X hl respectively, in analogy to the definition 
in Ref. l2ll . The angle between the spin quantization axis 
of the left lead, hl, and the spin quantization axis of the 
right lead, fipt, is given by (j). We show a sketch of the 
coordinate system and of the magnetization directions of 
the leads in Fig. [TJ As spin quantization axis of the dot 
we take the z-axis of this coordinate system. With this 
choice the tunneling Hamiltonian reads 



unncl 



^/2 



EK.+ (^^''X+^-^''X) (4) 



h.c. 



The tunnel matrix elements, Vi^it) and VR(t), can be 
both time dependent. The generalized tunnel rates are 
defined by F, (t, t') = \ 2W* (t) (f) p^,, = 

5 I]cr=± Ta.a (i, t') and Ta{t) = Ta{t,t). Here p^^a is 
the density of states of the spin species a in lead a, which 
is supposed to be constant. The spin polarization of lead 
a is defined as 



Pa = — 

Pa+ + Pa- 

and it can take values between and 1. 



(5) 



It consists of the Hamilton operators for the dot, for the 
left (L) and right (R) lead and for electron tunneling 
between dot and leads. The single-level quantum dot is 
described by the Hamiltonian 



dot 



(2) 



where the operator dl^{d^) creates (annihilates) an elec- 
tron with spin a =t, i on the dot and rio- = d'^^d^ is the 
number operator for electrons with spin a. The strength 
of the Coulomb interaction between electrons on the dot 
is denoted by U, which can be arbitrarily large. The en- 
ergy level e<j(0 = e<j + 5eo-(i) of the dot can vary in time. 
In the following, we assume the dot level to be spin de- 
generate, i.e. ef(t) = €i{t) — e{t). The Hamiltonian of 
the lead a, with a = L, R is given by 



Icad.f 



E ^'"'^l' 



k(T akcr • 



(3) 



ka 



We choose the spin quantization axis of lead a along 
the direction of its magnetization, ha. The spin a of an 



B. Real-time diagrammatic approach 

The Hilbert space of the single-level quantum dot 
is four dimensional, and it is spanned by the states 
X — ^t'I T It d (empty dot, singly occupied dot with spin 
up, singly occupied dot with spin down, doubly occupied 
dot). On the other hand, the non-interacting leads at- 
tached to the dot have a large number of degrees of free- 
dom and act as baths. Hence, we can trace them out to 
obtain an effective description of the quantum dot. The 
dot dynamics are fully described by its reduced density 
matrix, pdot, with matrix elements P^^^ = (x2|pdot|xi)- 
We introduce also the notation — for the diagonal 
matrix elements (probabilities). The time evolution of 
the reduced density matrix is governed by the general- 
ized master equation 



d 

dt' 



^Pil it) 



+ 



Xi ,X'2 



(6) 



dt'W''''',' (t,t')P^' (t') 
X2X2 ^ ' ' X2 ^ ' 
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The kernel W^^^,^ (^j connects the states x'l and X2 
time t' with the states xi a-nd X2 at time i. It is useful 
to define the vector of the average occupation probabil- 
ities P = (Po,Pi,Pd) - iPo,P^ + Pi,Pd) and the spin 



expectation value, in units of h, S 
components are given by 



{Sx, Sy, Sz), whose 



Pl+Pl 



; Sy 



(7) 



Since we consider spin-degenerate dot levels, i?| — Ei, 
we can drop the first term on the r.h.s. of Eq. ([5]). 

We are concerned with adiabatic pumping, where no 
transport voltage is applied across the dot. The leads are 
therefore described by the same Fermi function /(w). As 
a consequence, the instantaneous current through the dot 
vanishes, and we need to consider the first adiabatic cor- 
rection in order to obtain the pumping current through 
the dot. We perform an adiabatic expansion along the 
lines of Ref. [l9l We start by performing a Taylor ex- 
pansion around the time t of P(t') appearing inside the 
integral on the r.h.s. of the generalized master equation 



dt 



E 

x'l ,xi 



Pit (*) 



it' 



t) -p^/ 

' dt X2 



it) 



(8) 



The subscript t in Eqs. ^ and (fTO)) denotes the time 
with respect to which the adiabatic expansion is per- 
formed. This time t enters the respective quantities 
parametrically; both the instantaneous and the adiabatic 
correction to the kernel are functions of the time differ- 
ence {t — t'). At this stage, it is convenient to intro- 
duce the zero-frequency Laplace transform of the kernel 
lyXiXiV ' ^ rt rf^/ (^xiXiV ' ut^y In order 

to evaluate the kernel of the master equation we perform, 
on top of the adiabatic expansion, a perturbation expan- 
sion in the tunnel coupling F. In the following we take 
into account processes in first order in the tunnel cou- 
pling. This approach is valid in the weak-coupling limit, 
i.e. k^T ^ F. At the same time the condition for adia- 
baticity, F ^ fi, needs to be fulfilled. The instantaneous 
occupation probabilities and their adiabatic corrections 
obey the equations 



= 



V X2X2 / t 

x'i.x'2 



V x2X2yt \ X2J 



, (i,0) 



IfpX. 

dt ^ 



X2 )t 



Xi ■.X2 



X2X2 / 1 



(■',1) 



i^ii) 



{a, 



(11) 



(12) 



This expansion is justified by the fact that the response 
time of the system is much smaller than the timescale 
of the parameter variation in time. We then expand the 
kernel itself as 



W''''',' (t,t') 
X2X2 ^ ' ' 



t'). 

(9) 

The superscript i{a) denotes the instantaneous contribu- 
tion (its adiabatic correction). The instantaneous con- 
tribution corresponds to freezing all parameters to their 
values at time t, i.e. X{t) X{t). The adiabatic cor- 
rection is obtained by linearizing the time dependence of 
the parameters, i.e. X{t) X{t) + {T — t)d/dTX{T)\r=t, 
and retaining only first-order terms in the time deriva- 
tives. Finally, we perform the adiabatic expansion of the 
elements of the reduced density matrix 



P}2 it) 



(pXl 

\ X2 



(^X"2^ 



(10) 



The number in the superscripts designates the order in 
the perturbation expansion in the tunnel coupling. The 
fact that we find elements of the reduced density matrix 
in minus first order in the tunnel coupling is consistent 
with our perturbative scheme, as those terms are propor- 
tional to f2/F, which is small in the adiabatic limit. The 
evaluation of the matrix elements of the kernel is done 
using a real-time diagrammatic technique, which was de- 
veloped in Ref. [13, extended to systems containing fer- 
romagnetic leads in Ref. [2l|, and extended to adiabatic 
pumping in Ref. [T^ The adiabatic correction to the 
matrix elements of the kernel does not appear in lowest 
order in the tunnel coupling, as considered in this pa- 
per. The equations for the instantaneous probabilities 
and for the adiabatic correction to the probabilities can 
be summarized as 




l-/(e) 

- [1-/(6) +/(£ - 

fie + U) 



U)] 2 [1 -/(£+[/)] 

-2[l-/(e + C/)] 



(13) 




l-/(e) 
fie) -fie 
-fie + U) 



U)] ^2F„S.p„ 
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Similarly, the equations for the expectation value of the spin read 



4s = 

dt 



f{e)Po - ^ [1 - f{e) ~ /(e + U)] Pi - [1 - f{e + U)] Pa 



- r[l-/(e) + /(e + C/)]S + Sx^B„ 



(14) 



where we introduced the notation Pq. = pai^a- The 
interaction-induced exchange field or effective -B-field ap- 
pearing in Eq. (|14p is given by the principal- value integral 



(15) 



The instantaneous elements of the reduced density ma- 
trix are obtained by setting the left hand side of Eqs. 

and dH]) to zero and by assigning superscripts (z, 0) 
to the vectors P and S on the right hand side of the equa- 
tions. The first adiabatic corrections to the elements of 
the reduced density matrix are obtained by assigning to 
the vectors P and S the superscripts (i,0) on the left 
hand side of the equations and the superscripts (a, — 1) 
on the right hand side. 



III. RESULTS 

Starting from the master equation we compute both 
the instantaneous matrix elements of the reduced density 
matrix and their first adiabatic correction. These are 
needed as an input for calculating the spin and charge 
current. 



A. Dot occupation and spin 

The fact that no transport voltage is applied to the 
system has important consequences on the occupation 
probabilities and the expectation value of the spin. To 
lowest order in the tunnel- coupling strength F, the in- 
stantaneous occupation probabilities are given by their 
equilibrium values, i.e. by the Boltzmann factors of the 
respective states 



X 



z 



(16) 



where (3 — l/k-oT is the inverse temperature, E^{t) the 
energy of the dot state X: and Z the partition function. 
The spin expectation value vanishes,— 



(17) 



When considering only first-order tunneling processes, 
the adiabatic correction to the reduced density matrix is 
linear in fl/T. While the spin polarization of the leads 
has no influence on the instantaneous probabilities, the 



situation is different for the adiabatic correction, which 
reads 



p(a, 



with the 



.Q 



with 



dt 



charge 
1 



it) 



(18) 



relaxation time given by T^i{t) 
r{t) [1 + f{e{t)) - f{e{t) + U)] (the 



derivation of the expression for is given in Appendix 
lAl ) For vanishing polarization in the two leads this result 
coincides with that obtained for an N-dot-N syste m^^i'^^ 



p(a -1) ^ 



dpiifi) 1 



1 



dt r{t)l + fie{t))-f{e{t) + U) 



(19) 

We find non-vanishing contributions to the off-diagonal 
terms of the reduced density matrix, which vanish for 
zero polarization in the leads. The spin expectation value 
reads 



2 



(i,0) 



dt 



'rel 



it) 



K^o) 



where the spin relaxation time T^^iit) is given by 

{rrUt)y' = m [1 - fm) + fm + U)] (the derlva- 
tion of the expression for r^^j is given in Appendix [A|). 
Notice that the first adiabatic correction is the leading 
contribution to the expectation value of the dot spin. 
Furthermore, a time-dependent dot spin can be accumu- 
lated only by varying in time the occupation of the dot. 
The adiabatic correction to the spin component is parallel 
to the exchange field, introduced in Eq. (fT5| . Therefore 
no precession of the spin around this field takes place. 
This is different from the case of a time-independent but 
biased spin valve 

Finally, we remark that the limit where both leads are 
fully polarized along the same magnetization axis, (f) — 
and PL = PR = 1 is ill defined; in fact, in this case 
the life time of a minority spin in the dot diverges and 
consequently in order for the adiabatic expansion to hold 
the pumping frequency needs to be zero. 



B. Pumping Current 

The results for the dot occupation probabilities and the 
expectation value of the spin on the dot serve to calculate 
the pumping current. Using a similar approach as for the 
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generalized master equation we write the current into the 
left lead as 



where 



Xl,X2,Xi,X 



(^xf?) ('''^') "^'^ processes, describing 

transitions where the difference of the number of elec- 
trons entering and leaving the left lead is equal to the 
integer number q. 

We compute the first order adiabatic correction to 
the current including only first-order tunneling processes. 
We find 



(21) 

and 



Aa,0) 



(a,-l) 



Xl>X2,Xl:X2 



In the following we suppress the superscript (a, 0) for the 
current, since the instantaneous current is always zero 
and J^'^'"-' (t) is therefore the dominant contribution. The 
current is of zeroth order in the tunnel coupling and pro- 
portional to the pumping frequency f2. To this order 
in the tunnel- coupling strengths, the pumped current is 
non-vanishing only if the dot level position is one of the 

pumping parameters, since [P^'^j is proportional 

to the time derivative of the dot level position. We find 
for the pumping current 



~^ di Yitf 



dt 



(22) 



where we have defined the quantity 7r(i) = 

J2a^a{t)Pa/T{t) = Y.a'^a{t), which dcpCuds On 

time via Ta(t). The pumped current consists of two 
terms of different origin: the first one is independent 
of the lead polarizations, and can be interpreted as 
arising from a peristaltic mechanismji^ the second 
term depends on the polarizations and can be seen 
as arising from the relaxation of the accumulated 
spin on the dot. The latter contribution due to spin 
relaxation can be either positive or negative depending 
on the polarization strengths, polarization directions 
and tunnel coupling to the different leads. This means 
that the time-resolved current can be enhanced with 
respect to the non-magnetic case, in contrast with the 
spin-valve effect in a time-independent system, which 
always leads to a current suppression. Similarly, we will 
see later that also the charge pumped through the dot 
per period is always suppressed due to the polarization 
of the leads. Therefore, this inverse spin-valve effect is 
observable only in the time-resolved current response. 
The effect could be investigated experimentally by 



means of time-resolved measurements or by rectifying 
the current response. Reduction or enhancement of the 
current can be achieved by tuning the tunnel coupling 
or the lead magnetizations. 



C. N-dot-F: Spin Pumping 

We now turn our attention to spin pumping in a setup 
where only one of the leads, the right one for the sake 
of definiteness, is ferromagnetic. We calculate the spin 
pumped in the unpolarized left lead. 

The instantaneous contribution to the reduced density 
matrix is independent of the polarization and therefore 
remains unchanged. The first adiabatic correction for 
the occupation probabilities and the dot spin are given 
by Eq. ^ and Eq. respectively, with pl = 0. 

The charge current through such a N-dot-F system can 
be obtained directly from Eq. ([22]) by setting ph — and 
it reads 



(23) 



For calculating the spin current, we chose a global 
spin-quantization axis parallel to the magnetization of 
the right lead. In this basis the reduced density matrix 
does not have any off-diagonal terms. We find for the 
spin current 



hd{r 



rK{t)pMt) 



dt F(i)2-p|rR,(t)2 



(24) 



The ratio of the time-resolved spin and charge currents, 
Eq. (1221) and Eq. dSJ), reads 



h/2' 



(25) 



This ratio is, in general, time dependent. The time- 
resolved spin current is smaller than the time-resolved 
particle current at any time t. The ratio of these two 
currents is always positive, implying that spin and charge 
fiow in the same direction, as expected. The situation is 
different, for the spin (in units of h/2) and the charge (in 
units of — e) pumped per period. We calculate the pumped 
charge and spin for the following two choices of pumping 
parameters, {FL,e} or {FR,e}, in bilinear response, i.e. 
we calculate the pumped charge and spin per infinitesi- 
mal area in parameter space. The result for the ratio of 
the pumped spin (in units of h/2) per period, Ns, and 
the pumped charge (in units of — e) per period, N, reads 



Ns 
N 



~PR 



l+pj (fR/f)^-2(fR/f) 

i+pl{tn/r)'-2{tn/r)pl 



(26) 



It turns out that the efficiency of the spin pump does 
not depend on which pair of pumping parameters one 
chooses. In Fig. Ufa), we plot the ratio of pumped spin 
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FIG. 2: (a) Ratio of the pumped spin per period (in units 
of h/2) to the pumped charge per period (in units of — e) as 
a function of the relative tunnel-coupHng strength Fr/F, for 
different polarizations of the right lead, (b) Ratio between 
the linear dc spin conductance and the linear dc conductance 
as a function of the relative tunnel-coupling strength Fr/F, 
for different polarizations of the right lead. 



to pumped charge as a function of the relative tunnel- 
coupling strength r^/r, where the bar indicates time- 
averaged quantities, for different values of the polariza- 
tion of the right lead. The absolute value of the ratio is 
maximally equal to one in the case of full polarization of 
the right lead. For pR < 1, this ratio goes from — pR for 
vanishing Fr to pR for vanishing Fl changing its sign for 



Fr 
F 



(27) 



This is a very intriguing result, which implies that the 
respective direction in which spin and charge are pumped 
depends on the coupling to left and right lead. 

The average pumped charge and the average pumped 
spin can have opposite signs, while the time-resolved spin 
and charge currents flow in the same direction at any in- 
stant of time, due to the fact that the ratio of the time- 
resolved currents, Eq. P5|) . is itself time dependent. To 
elucidate this, in Fig. [31 we plot the time-resolved spin 
and charge currents as a function of time for a config- 
uration, where the pumped spin and charge per period 
have different signs. Note that the charge current has a 
positive average and the spin current has a negative av- 
erage, while both currents flow in the same direction at 
any time. 

We now contrast the results for the pumped spin and 
charge with the dc transport properties of the N-dot-F 
system. We find for the spin and charge currents 



1^ 



h/2 



[1-/L(e) + /L(e + (7)]FLPR 



[1 - /L(e) + /L(e + C/)]Fl + [1 - /R(e) + /R(e + U)]T^{1 ~ pi) 



r 



(28) 



which, in the linear response regime, yields for the ratio 
of the spin to the charge conductance 



G 



h/2 



Fl + Fr(1-p|) 



(29) 



The linear conductance ratio is shown in Fig. [^^b). Its 
behavior is completely different from that obtained by 
pumping. First, the spin polarization decreases as a func- 
tion of Fr/F and, second, it stays always positive. 

Finally, we consider the spin accumulated on the dot 
in one pumping period. We find two different results 
depending on whether {Fl, e} or {Fr, e} are the pumping 



parameters: 

(S("-^))t({Fl,4) = 

(S(''-i))t({Fr,6}) = 



r] d{n) 
7] d{n) 



^(1 



' rcl " 



(30) 

-PR , 

(31) 

where the area of the cycle in parameter space, 77, is de- 
fined as 77 = dt^STi^ and 77 = dt^STB. for the first 
and second equation respectively. In the case of pumping 
with {Fr, e}, i.e. when the coupling to the ferromagnetic 
lead is time dependent, the average spin changes sign at 
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FIG. 3; Time-resolved spin and charge currents as a function 
of time. The value of the parameters used for this plot are: 

PR = o.99,fL = Fr, e = -r, u = lor, |5rL|/f = \5e\/r = o.i 

and ksT = F. 



the same values of Fr/F at which the the ratio /l/^l 
changes its sign. On the contrary, when pumping with 
Fl and e, the average spin polarization of the dot does 
not change sign as a function of Fr/F, while the ratio 



stiU does. 



D. F-dot-F: Spin- Valve Effect 

We now consider the spin-valve setup with both 
leads having arbitrary spin polarizations. We com- 
pute the number of pumped charges per period, N — 
— i dtli^ (t) , in bilinear response in the pumping pa- 
rameters. For the pumping cycle defined by e(t) — 
e + 6e{t) and Ti^{t) = Ti^ + STi,{t), the number of pumped 
charges per period reads 



N 



V- 



d{n 



de 



(32) 



d ( rLEc.fa-rLPLEar„p„ 



f^-(EofaPa) 



with rj = dt^STi^ being the area of the cycle in pa- 
rameter space. Notice that the charge number in Eq. 
([32)1 is a product of two terms, where one contains ef- 
fects of interactions and another one effects of the leads' 
magnetization. 

In the following we show results for the case that both 
leads have the same spin polarization strength. This cor- 
responds to the experimentally relevant situation that 
both leads are realized with the same ferromagnetic ma- 
terial. In Fig. [31 we show the pumped charge as a func- 
tion of the level position for different values of the angle 
between the directions of the magnetizations of the two 
leads. The pumped charge shows a peak when the en- 
ergy e or e + U are close to the Fermi energy, similarly to 
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FIG. 4: Pumped charge as a function of the average level 
position, e, for different values of the angle between the mag- 
netizations. The polarizations in the leads are pl = PR = 0.8. 
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FIG. 5: Pumped charge as a function of the angle between 
the magnetizations of the leads, for different polarization 
strengths Pl = Pr = p. This result does not depend on the 
level position and the interaction strength. 



pumping through a quantum dot contacted to two non- 
magnetic leads (N-dot-N).^^ As far as the dependence on 
the angle between the magnetization of the two leads 
is concerned, for (/) € [0, tt], the charge is monotonically 
suppressed for increasing until a minimum is reached 
for (/) = TT as in the usual dc spin-valve effect. The full 
(/)-dependence of the pumped charge is shown in Fig. [5l 
where we plot N{(l))/N{(t) — 0). This result does not 
depend on the value of the level position and of the in- 
teraction strength, since the dependence on e and U can- 
cels out when we divide by N{(t) = 0). We notice that 
the suppression of charge pumping is stronger for higher 
lead polarizations. Furthermore, the more the lead po- 
larization is increasing the stronger the behavior of the 
pumped charge as a function of the angle deviates from 
a cosine law. 

In Fig. [6l we show the pumped charge as a function of 
the lead polarization strengths for different values of the 
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FIG. 6: Pumped charge as a function of the polarization 
strength p = pl ~ Pk, for different values of the angle be- 
tween the magnetizations. This result does not depend on 
the level position and the interaction strength. 



regime of weak tunnel coupling between dot and leads, by 
means of a real-time diagrammatic approach. In the case 
that only one lead is ferromagnetic, we have computed 
the spin injected in the non-magnetic lead by pumping. 
Wc have found that, depending on the relative strength 
of the tunnel coupling to the leads, spin and charge can 
be pumped, on average, in opposite directions. For the 
case when both leads are polarized, we have found a sup- 
pression of the pumped charge by means of the spin- valve 
effect and determined the average spin accumulated on 
the dot during one pumping cycle. 
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angle between the magnetization directions. This plot 
confirms that the pumped charge decreases for increasing 
spin polarization of the leads. The charge suppression is 
strongest when (j) is near to tt. Independently of the angle 
between the polarization axis of left and right lead, the 
pumped charge goes to zero for fully polarized leads. It is 
important to point out that this last property depends on 
the order in which limits are taken, since the two limits 
(j) —^ and PL — PR ^ 1 do not commute, as was already 
mentioned in section IIII Al In fact, comparing Fig. [5] 
and Fig. [S]we notice that in the first case the charge is 
maximal for (jj = even for the polarization increasing 
towards one, while in the second case for pl = PR = 1, 
the charge is maximally suppressed even for </) going to 
zero. 

The spin, which is accumulated on the dot during one 
pumping cycle is given by 



rj d{n 



de 



'rcl 



fpi 



r^-(EaP«ra) 



(33) 



(r^-(E.p«f„) 



where the pumping parameters are chosen to be Fl and 
e. The result for pumping with Fr and e is easily ob- 
tained by swapping the indices L and R. Depending on 
the spin polarization of the leads and on the values of 
tunnel-coupling strengths, the average spin on the dot 
can point along any direction in the plane containing the 
magnetizations of the leads. 



IV. CONCLUSIONS 

We have investigated adiabatic pumping through a 
single-level quantum dot with ferromagnetic leads in the 



APPENDIX A: RELAXATION TIMES 

In this appendix, we calculate the spin and charge re- 
laxation times. In order to calculate the spin relaxation 
time, we consider the case when the charge on the dot is 
in equilibrium and the occupation probabilities are there- 
fore given by the Boltzmann factors. Then Eq. (fT4)) 
simplifies to 



dS 



^F[l-/(e) + /(e + C/)]S 



(Al) 



where we also made use of the fact that the spin is always 
parallel to the exchange field. The spin relaxation time 
is therefore given by 



1 



1 



ri-/(e) + /(£ + [/) 



(A2) 



In order to calculate the charge relaxation time, we con- 
sider Eq. (jl3p . where we take the spin in equilibrium, 
such that S = 0. Then, we find for the dot occupation 
number 



d{r 



dt 



= F[2/(6)Po-(l-/(e)-/(e + i7))Pi 

-2{l-f{e + U))Pd] . (A3) 



Taking into account that the sum over the occupation 
probabilities has to be equal to one at any instant in 
time, we find 



d{n) 
dt 

where ( 

dot. The charge relaxation time is therefore given by 



/cq 



-F[l + /(£)-/(£ -I- f/)]((n)-(n),q) , (A4) 
is the equilibrium occupation number of the 



1 



'rel 



1 



Fl + /(e) -/(£ + [/) 



(A5) 
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Both relaxation times depend strongly on the position leads and on the strength of the Coulomb interaction, 
of the dot level with respect to the Fermi energy of the 
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This statement holds, also in the case of equally and fully 
polarized leads. The probability for a minority spin to enter 
the dot goes to zero, but also its probability to leave the 
dot does, once the minority spin is on the dot. 
Please note that the extra factor 1/2 in the formula in 
Ref. is a misprint. 



